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GLOBAL WEAK SOLUTIONS TO THE 
RELATIVISTIC VLASOV-MAXWELL 
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Abstract. In their seminal work [3], R- DiPerna and P.-L. Lions estab- 
lished the existence of global weak solutions to the Vlasov-Maxwell system. In 
the present notes we give a somewhat simplified proof of this result for the rel- 
ativistic version of this system, the main purpose being to make this important 
| result of kinetic theory more easily accessible to newcomers in the field. We 

show that the weak solutions preserve the total charge. 

f^) , Key words. Relativistic Vlasov-Maxwell system; global weak solutions; col- 

■ lisionless plasma. 
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1 Introduction 

. When a plasma is sufficiently rarefied and/or sufficiently hot like in the solar 

wind or in a powered-up fusion reactor collisions among the plasma particles are 
fH . sufficiently rare to be neglected. The only interaction among the particles then 

is through the electromagnetic fields which the particles create collectively. For 
the sake of simplicity we restrict ourselves to a plasma consisting of just one 
particle species, say, electrons, and we allow for the possibility that the particles 
move at relativistic speeds. The time evolution of the plasma is governed by 
" the relativistic Vlasov-Maxwell system: 

d t f+p-d x f + (E+pxB).d p f = 0, (1.1) 

dtE-cmlB = -4wj, d t B+cui\E = 0, (1.2) 
div E = 4irp, divB^O, (1.3) 

P(t,x)= f(t,x,p)dp, j(t,x)= pf(t,x,p)dp. (1.4) 



Here / = f(t,x,p) denotes the density of the particles on phase space, t Gl, x,p€ 
R 3 stand for time, position, and momentum, 

P 

p:=- 
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is the velocity of a particle with momentum p, E — E(t,x) and B — B(t,x) are 
the electromagnetic fields, p — p(t 7 x) and j=j(t,x) denote the spatial charge 
density and current, and units are chosen such that all physical constants such 
as the speed of light and the charge and rest mass of an individual particle are 
normalized to unity. The analysis can immediately be adapted to a plasma with 
several species of particles. 

We are interested in the Cauchy problem for the above system, i.e., in the 
existence of solutions satisfying the initial conditions 

o o o 

/|t=o=/i E\ t=Q = E, B\ t=Q = B (1.5) 

where the initial data satisfy the constraint part (|1.3|) of the Maxwell equations. 
Global existence and uniqueness of sufficiently smooth solutions to this initial 
value problem is an open problem. Local existence and uniqueness of classical 
solutions for smooth, compactly supported data was established in These 
solutions can be extended globally in time provided the momentum support can 
be controlled, which has been done for data which are small J3| or close to 
neutral 7 or close to spherically symmetric J22]. I n lower dimensions global 
classical solutions exist for general data |%1 131 1101 111) . Different approaches to 
the result in |12| were recently given in I16| . There is as yet no indication 
that classical solutions for general data in three dimensions develop singulari- 
ties. Nevertheless, it is natural to weaken the solution concept in order to obtain 
global solutions. This was done by R. DiPerna and P.-L. Lions [3J. The authors 
restricted themselves to the non-relativistic Vlasov-Maxwell system where p is 
replaced by p and pointed out that their arguments apply to the relativistic 
case as well. Many results on the Vlasov-Maxwell system are reviewed in the 
monograph jS], in particular, global weak solutions are discussed following [3J, 
but using also arguments introduced in |17j . The techniques in PJ are closely 
related to those used by the same authors in their seminal work on the Boltz- 
mann equation and later these techniques have been adapted to a variety 
of other problems in kinetic theory; we mention as being related to the 
present notes. 

Given the fact that the global existence and uniqueness problem for the 
Vlasov-Maxwell system in three dimensions is still open and the fact that in 
recent years many young researchers have joined the field it should be useful to 
give a simplified proof of the result of DiPerna and Lions. We emphasize that 
all the essential techniques we are going to use are introduced in [3] , but some 
non-trivial concepts and arguments from |3J are avoided. We concentrate on 
the relativistic version of the system, firstly, because being Lorentz invariant it 
seems better justified from a physics point of view, and secondly, because the 
non-relativistic case was considered in detail in |3J E] ■ In passing we note that 
there are global existence results for classical solutions of the Vlasov-Poisson 
system [^[20||23] but not yet for the so-called relativistic Vlasov-Poisson system 
where p is replaced by p in the Vlasov equation, cf. UJ. 

We now discuss how the paper and the proof proceed and how our ver- 
sion differs from the original one. In the next section we recall various a-priori 
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bounds resulting from conservation of energy and conservation of phase space 
volume by the characteristic flow of the Vlasov equation. Then we introduce 
a regularized version of the system which has global in time, smooth solutions. 
It will be important that these regularized solutions exist on the whole time 
axis. In |3] the system was regularized by adding a sufficiently large power of 
the Laplacian to the evolution part of the Maxwell equations l|1.2|l . This de- 
stroys the time reversibility of the system and technically impedes the straight 
forward application of the velocity averaging lemma discussed below. We pro- 
pose to regularize the system by smoothing the current j so that conceptually 
we remain closer to the unmodified system. This regularization is due to JS] 
and was used in the context of weak solutions in |SJ El D3| , but its technical 
advantages were not fully realized. Given a sequence of solutions to regularized 
Vlasov-Maxwell systems along which the regularization vanishes in the limit we 
show that the a-priori bounds hold uniformly. Hence we can extract a weakly 
convergent subsequence whose limit is the candidate for the desired weak solu- 
tion of the unmodified system. The main difficulty lies in passing to the limit 
in the nonlinear term in the Vlasov equation. This difficulty can be overcome 
since additional compactness of the approximating sequence is provided by the 
velocity averaging lemma. Its application is discussed in Section 5, and to make 
the present notes self-contained we give a proof of the relevant relativistic ver- 
sion in an appendix. Since the proof rests on Fourier transforming the Vlasov 
equation with respect to space and time it is advantageous that our approximat- 
ing solutions are defined on the whole time axis. In the set-up of [3] a certain 
cut-off and extension maneuver was necessary. A second problem lies in passing 
to the limit in moments of / like p and j, the difficulty being that the relevant 
weights in p are not test functions. As opposed to |Hj we derive weak conver- 
gence of p and j directly from a-priori bounds for these quantities. In all this a 
minimal requirement is that the initial data have finite energy and finite total 

o 

charge. The original proof assumed in addition that / is square integrable. We 

o 

assume that / is bounded, and we completely avoid the non-trivial concept of 
rcnormalization which was necessary in ^ — in the opinion of the present author 
the main motivation for studying weak solutions is not to allow for the greatest 
possible generality in the initial data but the lack of global existence results for 
stronger solution concepts. Having obtained a weak solution we examine some 
of its properties in Section 6. Sufficient regularity is established to make sense 
of saying that the solution satisfies the initial conditions. Then we prove that 
for almost all times t the total charge and more generally any L 9 -norm of f(t) 
equals its initial value, a result which was not obtained in [5] and for which 
we exploit the relativistic nature of the system. Since the Vlasov equation is 
a conservation law on phase space it is a desirable feature of any "reasonable" 
solution concept that solutions preserve the total charge. The energy at times 
t ^ is bounded by its initial value but is not known to be conserved, and neither 
are weak solutions known to be unique. 

To sum up, our proof is simplified compared to |3] in the sense that the only 
non-trivial tool employed is velocity averaging. The rest of the proof consists of 
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straight forward exploitation of straight forward a-priori bounds. 



2 Preliminaries and a-priori bounds 

We introduce the main conservation laws for the relativistic Vlasov-Maxwell 
system; the derivations are formal, and in which sense they hold depends on 
the type of solution under consideration. Writing the Vlasov equation in 
divergence form 

8 t f + div, [pf] + div p [(E(t, x)+pxB(t,x))f}=0 (2.1) 

and integrating with respect to p yields local conservation of charge, 

9 t p+div x j = 0. (2.2) 

On the level of the solutions of the characteristic system 

x~p, p~ E(t,x) +px B(t,x) 

of the Vlasov equation, conservation of charge is reflected in the induced flow 
on phase space being measure preserving and / being constant along the flow. 
A lengthy computation shows local conservation of energy, 

d t e + div x a = } (2.3) 

where the corresponding energy density and flux are defined by 

e(t,x) := j yflTWf(t,x,P)dp+^(\E(t,x)\ 2 + \B(t,x)\ 2 ), 

a(t,x) :— j pf(t,x,p)dp+—E(t,x)xB(t,x). 

J 47T 

These local conservation laws imply corresponding global conservation laws: 
Proposition 1 Consider a classical solution of the relativistic Vlasov-Maxwell 

o o o o 

system with compactly supported initial data /gC c (l ), f>0, E,B£C 2 ' 
satisfying the constraints $1.3p . As long as the solution exists its energy 

1 



y /T+Wf(t i x,p)dpdx+— J (\E(t,x)\ 2 + \B(t,x)\ 2 )dx 

as well as any L q -norm \\f(t)\\ q with 1 <q<oo are constant in time. 

We use the conserved quantities to derive bounds on p and j: For any R> 0, 



>+R^J Vi+Wfd P <c(^J VT+Wfdp^j ' , 



P(t,x) = / f(t,x,p)dp+ / f(t,x,p)dp 
l\p\<R J\p\>R 



-in 



< irR 6 \\f(t)\ 



4 



where for the last step we choose 



R= (y N/i+bF/dp) 



1/4 



and the constant depends on ||/||oo- Taking both sides of the estimate to the 
power 4/3, integrating in x, and using Proposition ^ we have the following 
a-priori bounds on p and hence also j which is dominated by p: 

Proposition 2 Along any solution as considered in Proposition^ 

l|p(*)IU/3,lli(*)ll4/3<C 

o 

where the constant C depends on the energy of the initial data and on ||/||oo- 

To conclude we observe that the constraints l|1.3l) propagate: If we have a 
solution of the system Ijl.lH , i|1.2D , i|1.4|) satisfying the constraints l|1.3|) initially 
then the constraints hold as long as the solution exists, since 

d t (div x E - Airp) = div x d t E - 4ird t p = div x (d t E + Airj ) = div x (curl^S) = 0, 

d t div x B = — divj, (curl^ E) = 0. 



3 The regularized system 

For a mollifier 

deC c °°(M 3 ), d>0, J d=l, deven (3.1) 

we consider the regularized relativistic Vlasov-Maxwell system (|1.1H , l|1.2|) , (|1.4|l 
where in the Maxwell equations (|1.2|l we replace j by d*j, the convolution re- 
ferring to x; recall that the constraints (|1.3|) propagate once they hold initially. 
Along a (local) solution of the regularized system charge is still conserved, hence 
/ \j{t) \ < / p{t) < C an d therefore all spatial derivatives of d* j are bounded uni- 
formly in t. This is sufficient to show that the regularized system has global- 
in-time solutions for initial data as specified in Proposition ^ details of the 
proof can be found in |15|. However, there is one difficulty with the regularized 
system: The time derivative of the energy is 

(j-E-(d*j)-E) dx 

which need not vanish. Hence, when constructing a sequence of approximating 
solutions by the regularization above we have to make sure we preserve the 
a-priori bounds, uniformly along the sequence. Let 

^(M^Lei^K 6 )!^, JJ ^T\pjg(x,p)dpdx<^ 
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be endowed with the weighted Z^-norm || • ||ki n with weight + \p\ 2 - We fix 
initial data 

/e^ in nL°°(R 6 ), E, B e L 2 (R 3 ) (3.2) 

satisfying the constraint (|f .3|l in the sense of distributions, and take sequences 
(/„)cC* c °°(R 6 ), bounded in L°°(R 6 ), and (E~), C C* C °°(R 3 ) such that for 
any pG [l,oo[, 

/„-►/ in4 in ni p (R 6 ), ini 2 (M 3 ). 
With d as above we let d n (x) :=n 3 d(nx), define 

E n := d n *E~^ E, B n :=d n *B~^B in L 2 (R 3 ), 
and denote by (f n ,E n ,B n ) the global solution of the regularized system with 

o o o 

the initial data (f n ,E n ,B n ), where j now is replaced by d n *d n *j n — the two 
d n 's are intentional — with j n defined in terms of /„ as in ljl.4|l . This solution 
exists for all £ € R by the reasoning above. By uniqueness, 

(E n ,B n )(t)=d n *(E~,B~)(t) 

where (E~,B~) solves the Maxwell equations l|1.2|l with initial data E~,B~ 
and current d n * j n . The energy 

JJ ^l + Wfn(t,x, P )d P dx+^ J (\E~(t,x)\ 2 + \B~(t,x)\ 2 )dx 

now is constant in time, its time derivative becoming 

J (jn ■ E n - (d n * j n )-E~)dx = 0. 

Since the modified energy defined above dominates the energy of the regularized 
solution we can use the arguments in the previous section to prove the following 
result, cf. Proposition [21 

Proposition 3 For initial data as specified in \3. £)) there exists a constant C > 
such that for any n£N the solution (f n ,E n ,B n ) of the regularized initial value 
problem described above exists for all i G R and satisfies the a-priori bounds 

Wfn WHkin, ||/»(*)||oo, \\E n (t)h, \\B n (t)\\ 2 , \\ Pn (t) || 4/3 , \\j n (t)\\ i/3 <C. 

4 The weak limit 

Proposition 4 There exist functions 

/ei 00 (R;ii in nL 00 (R 6 )), £,£e£°°(R;£ 2 (R 3 )), p, j G L°°(R ;j L 4/3 (R 3 )) 
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such that up to a subsequence 



/„-/mi 2 (Jxl 6 ), E n ,B n ^E,B inL^IxM 3 ), Pn ,j n ^ p ,j mL 4 / 3 (/xt 3 ) 
for any bounded interval JcR, 

/>0, p = / fdp, j= / pfdp a. e., 



and the Maxwell equations ifi.il)) . \ 1.31) as well as local conservation of charge 
hold in the sense of distributions. For almost all t G K the energy of 
(f(t),E(t),B(t)) is bounded by its initial value. 

Proof. The extraction of the weakly convergent subsequence is standard, and 
by a diagonal sequence argument this subsequence can be chosen independently 
of the interval /. Since /„ > for ji€N, the weak limit / is non-negative al- 
most everywhere. That the limits lie in the asserted function spaces is straight 
forward. For example let 4cl be bounded and measurable with Lebesgue 
measure X(A), and let R>0. Then 

// (7 VT+Wfdp+^(\E\" + \B\ 2 ))dxdt 
JaJ\x\<r\J\p\<r 87r ) 

f// (/ V^+Wfndp+^(\E n \ 2 + \B n \ 2 ) ) dxdt 
Jaj\x\<r \J\p\<r 87r / 



< liminf 



1 



<HA){\\f\U + —{\\E\\ 2 2 + \\B 
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Since i?>0 and Act are arbitrary the assertion on the energy follows, which 
implies part of the assertion on the weak limits f,E,B. 

For the moments p and j we can argue similarly, but we have to make sure 
that the weak limit of p n is the density p induced by / and analogously for j. 
To this end, consider a test function t|ieC c °°(Rxl 3 ) and some R>0. Then 

pfdp — j I ipdxdt 



pfdp — j ipdxdt+ pf dpipdxdt 

\p\<R I JJJ\p\>R 



lim I I \ I pfndp-jn i>dxdt+ / / / pfdpipdxdt 

\p\<R I JJJ\p\>R. 



n — >oo 



lim /// p{f - f n )dpipdxdt. 
l\p\>R 



The modulus of the latter integral can be estimated by 

M«i // [ V^+W(fn+f)dpdxdt<^- 



via the uniform bound on the kinetic energy, and since R > and the test func- 
tion ip are arbitrary the assertion for j follows. That local conservation of charge 
l|2.2[l and the second of the Maxwell equations (|1.2[1 hold in the sense of distri- 
butions is obvious. As to the first of the Maxwell evolution equations let tjj be 
as above. Then with the abbreviation 5 n :=d n *d n , 



5 n *j n ipdxdt- J J jtpdxdt = J J (5 n * j n - j n )tpdxdt + J J (j n - j)if)dxdt. 
The second term converges to zero, and the first can be estimated by 

IUn||iV3( su p pV ,)||^„*V J -V J ||4 

which converges to zero as well. Hence S n *j n — > j in the sense of distributions. 
The computation in Section 2 which showed that the constraints ()1.3|l propagate 
can be performed in the weak sense, and the proof is complete. □ 



5 Compactness via momentum averaging 

In the present section we want to show that the weak limit obtained in the 
previous one satisfies the Vlasov equation in the sense of distributions. Due to 
the non-linearity in the latter equation this is the crucial problem in the whole 
proof. We need to show that for any ^C c °°(lxR 6 ), 

(E n +px B n ) -dp4>fndpdxdt~^ JJJ (E+px B) ■ d P 4> f dpdxdt, 

possibly after extracting a further subsequence. By a well known density ar- 
gument it is sufficient to consider test functions which factorize: (j}(t,x,p) = 
(j)i(t,x)(j)2(p) with test functions 0i,02- The desired convergence will follow 
once we can show that for any bounded open interval Jcl, any S > 0, and any 
V' I =C^ (]R 3 ) up to a subsequence, 

J ^{p)fn(-,-,p)dp^ J ip(p)f(;-,p)dp strongly in L 2 (IxB s ), (5.1) 

where B$ C M 3 denotes the ball of radius S about the origin. That this is 
true is the consequence of the velocity-averaging lemma, a version of which was 
established in 0, cf. also ^1], and which for obvious reasons we prefer to call 
momentum-averaging lemma. We state it here and give a proof in an appendix: 




Lemma 1 Let R>0 and i\> € C^{Bji). There exists a constant C > such that 
for any functions h, go,gi € L 2 (R x R 3 x Br) which satisfy the inhomogeneous 
transport equation 

d t h+p-d x h = g + div p gi (5.2) 
in the sense of distributions we have 

J V>(p)Mv,p)dpeff 1/4 (RxR 3 ) 
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with 



<C(II^IU 2 (kxe 3 xb r ) + ||3olU 2 (---) + HgilU 2 (- 

#1/4 



Here iJ 1 / 4 denotes the usual fractional order Sobolev space denned in terms of 
the Fourier transform. In the proof one needs to Fourier transform the transport 
equation with respect to t and x so that it is essential that the equation holds 
for all t G R. Since we apply the lemma to the approximating solutions we avoid 
considerable technical complications by having the latter defined on the whole 
real line. 

In order to prove (|5.1(l we take an arbitrary open, bounded interval Id 
R, ■0GC' c oo (R 3 ), and R>0 such that supp-0 C Br. Moreover, we choose some 
test function CgC c °°(R) such that 0<C<1 and C|/ = l- Then for f n (t,x,p):= 
C(t)fn(t,x,p) the Vlasov equation in divergence form (|2.1|l implies that 

dtfn+P-d x f n = go + div p gi 

where g$ := ('f n and gi := (E n +px B n )(f n . By Proposition |3 the L 2 -norms of 
the latter functions over the domain lxl 3 x Br are bounded, uniformly in n. 
Hence by Lemma ^ we conclude that the momentum averages of the functions 
f n , formed with the test function tp, lie in the Sobolev space 7J 1 / 4 (RxR 3 ), 
and their corresponding iJ^-norms are bounded, uniformly in n. Over the 
open, bounded set I x Bs with S > arbitrary this Sobolev space is compactly 
embedded in L 2 , by choice of £ the functions /„ coincide with /„ on / x Bs , 
and hence along a suitable subsequence Eqn. I|5.1|l is established. 

We conclude that the weak limit (f,E,B) obtained in the previous section 
satisfies the complete relativistic Vlasov-Maxwell system on R x R 3 x '. 
sense of distributions. 



6 Continuity properties and conservation laws 

In order to be justified in saying that the weak solution (f,E,B) satisfies the 
initial conditions we need to establish some minimal continuity in t. We consider 
this problem for /, the arguments for E and B being very similar. We express 
dtfn via the Vlasov equation in divergence form, integrate in time, multiply the 
result by a test function 0eC^°(R 6 ) and integrate in x and p to obtain, after 
an integration by parts: 



f n (t,x,p)<j)(x,p)dpdx = J J f n (x,p)4>(x,p)dpdx 

\p-d x (f>+(E n +px B n )-d p (t)]f n (s,x,p)dpdxds. 
For t G R we define a distribution f(t) G X>'(R 6 ) by 
(f(t),4>) : = f(x,p)(j)(x 1 p)dpdx 
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ai 



t 

[p-d x (j)+{E+px B) ■ d P 4>] f (s,x,p) dpdxds. 



The mapping 1 1— > /(t) G 2?'(R 6 ) is continuous in the sense of distributions, /(0) = 

o 

/, and one can show that this is a representative of /. By a density argument 
1 1— > f(t) is also continuous with respect to the weak topology of L q (R 6 ) for any 
1 < q < oo. Since for any bounded measurable set Mcl 6 and s > the negative 
order Sobolev space W~ S ' 2 (M) is compactly embedded in L 2 (M) continuity 
with respect to the negative order Sobolev norm follows. 

We already noted that the total energy can at least not increase, and that 
local conservation of charge (|2.2(l holds in the sense of distributions. It is simple 
to conclude that for almost all t, / p(t) equals some fixed constant, but we assert 
that this constant really is the initial charge. Indeed, more is true: 

Proposition 5 The weak solution obtained above preserves all L q -norms of 
f(t), more precisely, for every l<q<oo, 

\\f(t)h = \\f\\ q for a. a. teR. 

Moreover, ti— >/(t) is a. e. strongly L q continuous in the following sense: There 
exists a set of continuity CcR such that the Lebesgue measure A(R\C)=0, 
OeC, and for any g€]l,oo[ the mapping 1 1— > f(t) S Z 9 (R ) is strongly con- 
tinuous. 

Proof. Since dtp n + div j„ = and \j n \<p n we have for every R > and t > 0, 



/ p n (t)dx = - p n (t)dS x + j 

J\x\>R+t J\x\=R+t J\s 



d t pn{t)dx 

\ = R+t J\x\>R+t 



d 

dt_ r 

p n {t)dS x - I divj n (t)dx 

x\=R+t J\x\>R+t 

(Pn(t) + V3n(t))dS x <0 

x\=R+t 

where v is the outer unit normal of the domain > R + t}. The analogous 
argument works for t < and the domain {|x| > R — t}. Hence 



/ p n (t)dx< / p n dx, teR, R>0, neN, 

J\x\>R+\t\ J\x\>R 



(6.1) 



|x|>-R+|t| 

O 

where p:= J f dp. We claim that 

J p{t)dx = J pdx for a. a. teR. (6.2) 
Let e > be arbitrary. Since p is integrable, we can choose R> such that 

pdx < e. 

\x\>R. 
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By the convergence of the initial data and l|6.1fl we conclude that 
/ p n dx< p n dx<e 

J\x\>R+\t\ J\x\>R 

for all i£R and all sufficiently large nSN. Let icM be measurable and 
bounded. Then 



pdxdt > 



pdxdt= lim 



AJ \x\<R+\t\ 



p n dx dt 



A J \ x\ <R+\t\ 



= lim 



pndx- 



Pndx 



dt>X(A) 



pdx — e 



<\x\>R+\t\ 

and for sufficiently large S > we have by monotone convergence 



pdxdt < 



I AJ\x\<S 

= lim 



pdxdt + X(A)e 

p n dxdt + \(A)e<\(A)[ I pdx + e 



iAJ\x\<S 

This implies that there exists a set M e cK of measure zero such that 

J pdx-e< J p(t)dx< J pdx + e, teR\M e . 
Hence H6.2JI holds on R\UfceN-Mi/fc. For general q£ [l,oo[ we define 
p q (t,x):= / f q (t,x,p)dp, j g (t,x):= pf q (t,x,p)dp, 



with analogous definitions for p q . n and j q , n - Since / and f n are bounded, these 
new densities converge in the same sense as the ones for q = 1 and the weak 
limits are indeed the g-densities induced by /. Moreover, dtp q . n + divj 9i „ = 
classically for every n, and \j q . n \ < Pq,n- Hence exactly the same argument as in 
the case q = 1 shows that 



p q {t)dx= I p q dx 



(6.3) 



for almost all t € R. If we pick a sequence qk — > oo we can choose the exceptional 
set of measure zero where does not hold uniformly in fc, and passing to the 
limit shows that the L°°-norm of f(t) is preserved in the same sense. 

Combining this with the weak continuity of the map tt—>f(t) or rather a 
suitable representative, the Radon-Riesz Theorem implies that for every q£ 
]l,oo[ one can choose a set of continuity C q with the asserted properties, and 
via two sequences q^ — > oo and qZ \, 1 we obtain a set of continuity C which 
works for all gs]l,oo[ simultaneously. □ 

Note that the relativistic nature of the system was exploited in the proof 
above. We are not aware of a proof of (|6.2|l in the non-relativistic case. 
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7 The result — statement and comments 



Collecting the results of the previous sections we arrive at the following existence 
result for global weak solutions to the relativistic Vlasov-Maxwell system: 

Theorem. For initial data /ei^ n nL°°(R 6 ), E, BeL 2 (R 3 ) which satisfy the 
constraints in the sense of distributions there exists a global weak solution 

of the relativistic Vlasov-Maxwell system, i.e., there exist functions 

/ei 00 (R;i ] ^ in nL 00 (R 6 )), £,5eL°°(R;L 2 (R 3 )), p, j G L°°(R ;j L 4/3 (R 3 )) 

such that (f,E,B) satisfy \1. 1)) — ]1.S\) in the sense of distributions with p, j 
defined in terms of f by j l-4\j . The function f is a. e. non-negative. 

The mappings tt—>f(t),E(t),B(t) are (after choosing suitable representa- 
tives) continuous with respect to the following topologies: the standard topology 
in the space of distributions £>'(R 6 ) or 2?'(R 3 ) respectively, the weak topology 
of L 2 , and the strong topology of W~ S,2 (M) for any s>0 and any bounded 
measurable subset M of R 6 or R 3 respectively. The mapping ti— » f(t) is a. e. 
strongly continuous into any L 9 (R 6 ), l<q<oo, in the sense of Proposition^ 
The initial conditions J^1.5}) are satisfied. 

At a. e. time t the total energy 

JJ^TTWf(t,x,p)dpdx+^- J (\E(t,x)\ 2 + \B(t,x)\ 2 )dx 
is bounded by its value att = 0. The total charge is conserved, 

J J f{t)dpdx = J J fdpdx for a. a. ieR, 
and the same is true for any L q -norm of f(t), l<q< oo. 

The following deficiencies of weak solutions are obvious: Firstly, uniqueness is 
not known. Secondly it is not clear whether energy is conserved and whether 
L^-norms of / are conserved everywhere. Moreover, their continuity with respect 
to t holds either everywhere in a rather weak sense or only a. e. in the strong 
L q sense. A less obvious disadvantage is the following: It is not clear that any 
weak solution in the sense of the theorem can be obtained via the particular 
regularization which we employed here, and it is conceivable that results for 
weak solutions depend on the way in which these are constructed. 



Appendix: Momentum averaging 

To prove Lemma ^ we let -0gC^°(R 3 ) with supp-0C-B_R and R>0. By the 
usual definition of the fractional order Sobolev spaces via Fourier transforms, 



H-,-,p)^(p)dp 



< 



h\\i + 



\I(r,0\ 2 (\r\ 1/2 + \e /2 )^dr (7.1) 
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where 



J(t,0:= / h(r,^p)^(p)dp 



(7.2) 



and h denotes the Fourier transform of h with respect to (t,x). Confusion with 
the notation p seems unlikely. By assumption, h satisfies the transport equation 
H5.2fl onRxl 6 in the sense of distributions, and hence 



i(T+p-£)h = g + div p gi. 



(7.3) 



This identity is only useful where the factor on the left hand side is away from 
zero. Hence we let (eC™(K) be such that 

0<C<1, sup P Cc[-2,2], C|[-i,i]=l 
and we split the integral 1)7. 2f> into the two parts 

'r+p-C 



HT,^,p)ip(p) 



\ K 



dp, 



1-C 



T+p-£ 



dp. 



where k > will be chosen appropriately in dependence of £. 
For almost all (r^)eR 4 , 



\h(T,0\<C\\h(T,£,- 



1/2 



1|t|<t|5|+2k 



(7.4) 



where r :— Rj \/l + R 2 — as in the rest of the argument, the L 2 -norm with respect 
to p refers to the ball Br with R > fixed, and 1 with some subscript is the 
"indicator function" of the set or condition in the subscript. To see ()7.4|) we 
apply the Cauchy-Schwarz inequality and observe that 



2 ( T+p-£, 



dp< I 1[- T -2 K -T+2 K ](pi\€\)dp; 



without loss of generality we may assume that (|£|,0,0). The integral on the 
right hand side vanishes if |r|>r|^| + 2K, and via a change of variables it can 
easily be estimated against which proves l|7.4|l . 

For the estimate of I2 the basic idea is to use l|7.3|l and to integrate by parts 
in the term containing the divergence of gi . With the abbreviation 



x(t,£,p) :=-0(p) 



i-C 



r+p-i 



1 



T+p-i 



it follows that 
h(r,0 = - 



9o(T,^,p)x(T,^,p)dp- / g 1 (T,£,p)-d P x(T,£,p)dp 



(7.5) 
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where this last identity holds a. e. onlxR 3 if k= 1 and a. e. on K x {|£| > 1} 
if k = |£| 1//2 . The first choice will be used for small |£| and the second one for 
large |£|. For any multi-index aeN 3 with length \a\ < 1, 



Hence 

H^xC-,e,-)lU<c(i+M) |^ (/^'i^i^^^) 7 • 

The remaining integral is straightforwardly estimated by C/k. In addition, for 
1 7" I > r \£, \ + 2k the origin does not lie in the domain of integration of this integral, 
and we can estimate it by integrating 1/w 2 over the domain. Hence, 

/ |m / i i \ 1/2 

\\DpX{T,t-)h<C{l + — j I l| T |< rK | +a<t — + l| T | >r | f | +aw -^— J . 

Let us abbreviate 

7V(r,e):=||Mr,^.)||2 + ||go(r,e,-)||2 + ||.9i(r,e,-)ll2. 

Combining l|7.4|l with l|7.5|l and the estimate for D^x we finally obtain the 
following estimate for J, 



\I(t,0\ < CN(t,0 



l|r|<r|g|+2n 



lei 



+ 1 kl>^l?l+2«; I i + — 



f) 2 








2 


i 




r 2 |£| 2 



1/2 



which holds a. e. on K x R 3 if k = 1 and a. e. on M x {|£|>1} if k=\£\ 1/2 - Hence 
f° r |£| < 1 wc take k = 1 to obtain 

1 1 X 1/2 



\I(t,£)\<CN{t,S) l|r|<r|«|+2^ + l|r|>r|?|+2 



and for |£| > 1 we take k= l^l 1 / 2 to obtain 

/ 1 |£| \ 1/2 

|i"(r,C)|<^iV(T,0 ^l|r|<r|5|+2|5|i/2|^p7J + l|r|>r|C|+2|C| 1 / 2 r 2_ ? ,2|^|2 J ' 

Now we split the right hand side of l|7.1|l as follows: 

IJ\I(r,mH\r\ 1/2 + \e /2 )drd^ J J \I\ 2 \r\^ 2 l m>1 drd^ 

+ J f \I\ 2 \T\ 1 ' 2 l^< 1 l\ T>r+2 dTd^+ ff llflT^l^l^r+zdrdt 

+||i/i 2 ier /2 i k |< 1 dr^+/y'i/i 2 ier /2 iici>idr^ 

=:Ai+A 2 +A 3 +Ai+A 5 . 
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Using the appropriate part of the estimate for I we find that 

A X +A 2 +A 5 <C j j N 2 drd^=C{\\h\\t + \\h\\lHh\\l) 

as desired, while the terms A3 and A4 can be estimated directly by ||/i(V,£,-)||2 
via the Cauchy-Schwarz inequality, and the proof of Lemma ^ is complete. 

Acknowledgment: The present notes originate in a seminar on kinetic 
theory at the University of Vienna. I would like to thank the participants of 
that seminar as well as Prof. R. T. Glassey for their interest in the simplifications 
I report here. 
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